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Abstract 

In this paper are investigated strictly self-similar fractals that are composed of an 
infinite number of regular star-polygons, also known as Sierpinski n-gons, n-flakes or 
polyflakes. Construction scheme for Sierpinsky n-gon and n-flake is presented where 
the dimensions of the Sierpinsky oo-gon and the oo-flake are computed to be 1 and 2, 
respectively. These fractals are put in a general context and Iterated Function Systems 
are applied for the visualisation of the geometric iterations of the initial polygons, as well 
as the visualisation of sets of points that lie on the attractors of the IFS generated by 
random walks. Moreover, it is shown that well known fractals represent isolated cases of 
the presented generalisation. The IFS programming code is given, hence it can be used 
for further investigations. 
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1 Introduction - the Cantor set and regular star-polygonal 
attractors 

A classic example of a strictly self-similar fractal that can be constructed by Iterated Function 
System is the Cantor Set [1]. Let us have the interval E = [—1,1] and the contracting maps 
Si,S 2 : E ^ E, = x/3 - 2/3, S 2 = x/3 + 2/3, where x e E. Also : S{S'^-\E)) = 
S^(E), S^{E) = E, where S{E) = Si(E)[J S 2 {E). Thus if we iterate the map S infinitely many 
times this will result in the well known Cantor Set; see figure This iteration procedure can 
be generalised by the following theorem [2] : 

Theorem 1. If we have Si,S n ■ \Si{x) — Si{y)\ < Ci\x — y\, Ci < 1, then 3 unique non-empty 
set E : E = U^iSi{E), hence invariant for the map S and E = n'^^iS’^{E) 
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Figure 1: Sketch of the repeated actions of the maps Si and S 2 on the interval E that result 
in the Canator Set, where the left and right arrows represent Si and S 2 , respectively On the 
right hand side of the intervals, the corresponding iterations of the map S and the number of 
the intervals of the particular iteration are shown. 

Using a polygon as an initial set for fractal generation is a well known technique since the 
most famous strictly self-similar fractal examples, the Cantor set and the Sierpinski triangle, 
consist of infinitely many line-segments and triangles, respectively. In the present paper the 
number of the vertices of the polygon will be increased to an arbitrary number n e N, n > 2. 
Thus, the fractal can consist of infinitely many pentagrams, hexagons, etc. Furthermore, the 
building regular polygons will be all {n/m} star-polygons [3] where n ^ 2 and m < n/2, 
n G N, and m G N. For our purpose we will take the unit circle and it will be divided in 
n equal segments. For example, in the case of pentagram, we have m = 2 which gives us 
{5/2}-polygon. Once we choose an {n/m}-polygon it can be scaled by a factor of P G (0,1) 
with respect to all of the vertices of the polygon. This will produce n new polygons similar to 
the initial one but scaled down by factor of P. Now, if we repeat the procedure for each one of 
these new polygons, another polygons will be created that will be P^ times smaller than the 
initial {n/m}-polygon inscribed in the unit circle. If the P is chosen carefully, after infinitely 
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many contractions, the result will be a strictly self-similar fractal composed by non-intersecting 
polygons. Thus, at the z-th contraction the defined n contracting maps are applied n* times, 
and by theorem [T] when z oo we will define infinitely many points of attraction. These points 
of attraction specify the attractor of the iteration procedure. This polygonal attractor is a 
fractal produced by the infinite contractions {n\ where i —>■ oo) of the initial polygon and it is 
self similar, i.e. it is composed of infinitely many polygons similar to the initial one. 

The present study is focused on non-self-intersecting fractals where the scaled copies of the 
initial polygons osculate with each other. This restricts the possibilities for P when Sierpinsky 
n-gon |1] for arbitrary n is constructed and a formula about P is derived. This important ratio 
is reported in two other places [1] and j5] where in the latter the proof is omitted. Moreover, 
in both manuscripts the authors did not prove that the Hausdorff dimension of the non-self- 
intersecting cxo-gon is 1. In the present paper, an original derivation of the equation for the 
scaling ratio P is presented. It is done in a very detailed way by using simple geometric laws 
which makes the result affordable even for high-school students. Also, the Hausdorff dimension 
of the oo-gon is shown to be 1. Furthermore, universal constructions for n-flakes are proposed 
for the cases when n is even or odd, and the Hausdorff dimension of the oo-ffake is proved to 
be 2. To this end, formulas for the scaling ratio and the rotation of the central polygon of the 
n-ffake are derived, which to the knowledge of the author have not been reported previously. 
Finally, it is shown that different initial polygons may result in an identical attractor when an 
IFS iteration scheme is applied and it is shown which are the main parameters that define the 
shape of the Sierpinski n-gon and the polyflake attractors. 

The paper is constructed as follows: in section]^ the parameters P and m are introduced and 
an important equation for the ratio P = P{n, m) is derived. In section a condition for m is 
obtained that ensures no self-intersection of the studied class of fractals. Then, two techniques 
for imaging IFS attractors are introduced and several Sierpinsky n-gons are computed together 
with their dimensions. In section]^ the possibility for additional scaling map that scales down 
towards the centre of the polygon is taken into account. Different constructions for the cases 
when n is odd or even are proposed so the resulted n-flake to be non-self-intersecting for 
arbitrary n. Also, a few interesting examples are given and the Hausdorff dimension of the 
(X)-ffake is computed. In section some well known fractals are shown to be a special case of 
the fractal generation scheme shown here. It is also explained why identical attractors may 
originate from different star-polygons and how we can exploit this feature. 
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Figure 2: Sketches of {n/3} and {n/m} star-polygons circumscribed in Sg, and (Sb respectively. 


2 The parameters P and m 

An important result of the present paper will be explained in this section. Here the scaling 
parameter P will be deduced from n and m. Therefore, P = P{n, m) is a specific scaling 
factor for the chosen initial {n/m}-polygon, where P does not depend on the diameter of the 
circumscribed circle. 

In figure]^ a sketch of a {n/3} star-polygon is shown where m = 3. Here, the vertices Aj for 
i = 1,..., n are the vertices of the {n/3} star-polygon and Oa is the centre of the circumscribed 
circle Sg, M is the intersection point of the secants A 1 A 4 and A^A^, H is the orthogonal 
projection of Oa on A 1 A 4 and L is the orthogonal projection of Oa on A 3 A 4 . Our purpose will 

be to find the ratio P = ^ ^ because MA 4 will be a line segment of the star-polygon resulted 

after the scaling of the initial polygon with respect to the point A 4 by factor of P. 

The initial polygon is circumscribed by Sg and Oa is its centre thus AA^OaA 4 ^ = 27r/n 

ZLOaA 4 = 7 r/n because A^Oa = A 4 ^ 0 a- Also, AAiOaA^ = dyr/n ZA 1 A 4 A 3 = 27r/n and 

AAiOaH = 37 r/n, since Sg with centre at Oa circumscribes Ai, A 3 and A 4 . Now we can deduce 

A 1 A 4 and A 4 M by the radius r of Sg] r = OgAi for i = 1 , ...,n : 

A 1 A 4 = 2 rsin( 37 r/n). In order to deduce A 4 M we will first find A 4 L. Thus, A 4 L = 

■ / I \ A ^ 4 ^ r sm(IT/n) . , , r . , , 

rsm TT/n) A 4 M = -;-- =-;-. Now we can substitute the values tor A 1 A 4 and 

^ ^ ^ cos(27r/n) cos(27r/n) 

,, . , ^ „ MA 4 sin(7r/n) 

MA 4 m order to find P — — 


A 1 A 4 2 cos(27r/n) sin(37r/n) 

We have just found P{n,3), now, let us do the same computations for any 1 < m < n/ 2 . In 
figure [^b) we can find the sketch of a {n/m} star-polygon where the vertices Aj for z = 1,..., n 
are the vertices of the star-polygon and Oi, is the centre of the circumscribed circle <Sb with 
radius r, M is the intersection point of the secants AiA^+i and AmA 2 m, H is the orthogonal 
projection of Ob on AiA^+i and L is the orthogonal projection of Ob on A^A^+i. Our purpose 

will be to find the ratio P = ——Thus, ALObAm+i = tt/u because A^O;, = A^+iO^, 

'‘^1 ''^ 777/-|-1 

AAiObAja = (2m —2)7r/n ZAiA^+iAm = (m— l)7r/n and AAiObH = mir/n, since «Sb with 
centre at Ob circumscribes Ai, A^ and A^+i- Now we can deduce AiA^+i and A^+iM by the 
radius r of (Sb: AiA^+i = 2r sin(m7r/n). In order to deduce A^+iM we will first find A^+iT. 
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Figure 3: IFS generated fractal sets made of points that lie on {9/2} and {9/3} star-polygons 
in subpanels (a) and (b), respectively. 


Thus, Am+iL = rsin(7r/n) ^ Am+iM = 
substitute the values for AiA^+i and MA 


Am+lL _ 

cos((m — l)7r/n) 
m+i in order to find: 


rsin(7r/n) 

- — - r—Now we can 

cos((m — IjTr/n) 


MAm+i _ _ sin(7r/n) _ 

AiAm+i 2 cos((m — l)7r/n) sin(m7r/n) 


( 1 ) 


3 Generation of fractals by using IFS 

The values for P obtained in the previous section will be used here for the computation of 
self-similar fractals that are IFS attractors. These attractors will be derived by the random 
walk/orbit method or so called chaos game [HI E]- First, we define a matrix that specifies how 
many points along the unit circle will be taken into account (n), and what is the contraction 
P(n, m) towards those points. For example, the matrices for the fractals in figures [^a) and 
l^b) are bellow; see Table 

This matrix is then plugged into the random generator, where the number of points that we 
want to map over the IFS attractor are specified; see the Appendix section for the MATLAB 
code. 


3.1 Condition for non-self-intersection 

In figure 1^ two examples of star-polygon fractals with initial {9/2}- and {9/3}-polygons clarify 

why the parameter m in the ratio ——is important when non-self-intersecting fractals are 

AiAm+i 

desired. We would like the self-intersection of the resulted sets to be prevented, thus, we will 
state the following theorem. 
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Table 1: The two matrices M{ 9 / 2 } and Mjg/s} where each of them defines nine contracting maps 
(every two rows define a map), needed for the random IPS procedure, resulting in the {9/2} 
and {9/3} fractal attractors shown in figure |^a) and|^b). 




^{ 9/31 


0.2831 

0 

0.8660254 

0.2578 

0 

0.8660254 

0 

0.2831 

-0.5 

0 

0.2578 

-0.5 

0.2831 

0 

0.9848078 

0.2578 

0 

0.9848078 

0 

0.2831 

0.1736482 

0 

0.2578 

0.1736482 

0.2831 

0 

0.6427876 

0.2578 

0 

0.6427876 

0 

0.2831 

0.7660444 

0 

0.2578 

0.7660444 

0.2831 

0 

0 

0.2578 

0 

0 

0 

0.2831 

1 

0 

0.2578 

1 

0.2831 

0 

-0.642788 

0.2578 

0 

-0.642788 

0 

0.2831 

0.7660444 

0 

0.2578 

0.7660444 

0.2831 

0 

-0.984808 

0.2578 

0 

-0.984808 

0 

0.2831 

0.1736482 

0 

0.2578 

0.1736482 

0.2831 

0 

-0.866025 

0.2578 

0 

-0.866025 

0 

0.2831 

-0.5 

0 

0.2578 

-0.5 

0.2831 

0 

-0.34202 

0.2578 

0 

-0.34202 

0 

0.2831 

-0.939693 

0 

0.2578 

-0.939693 

0.2831 

0 

0.3420201 

0.2578 

0 

0.3420201 

0 

0.2831 

-0.939693 

0 

0.2578 

-0.939693 


Theorem 2. If we have a strictly self-similar fractal set obtained as an attractor of IFS, where 
n-attracting points lie on so that they are the vertices of a {n/m} star-polygon, and where 
the attraction towards these points is P = P{n, m) given by equation 0 , then this fractal set 
is not-self-intersecting if and only if m E [n/4,n/4 + 1], which uniquely defines P for a given 
n. 

Proof. For clarity, one must look at figure]^ where with red is denoted the scaled down polygon 
towards A^+i, self-similar to the original one. Although, it has 9 vertices, it must be considered 
as n/m star-polygon, because we will only use geometrical properties that are independent of n 
and m. For this purpose we must find out the following angles: ZA'MA^+i and AObMA^+i- 
We already found that ALAm+iM = {m — l)7r/n and therefore 
AMAm+iOb = lObAm+iA2m+i = 7r/2 - Tx/n - {m- l)7r/n ^ 

^ ALAm+iA2m+i — 2(7r/2 — 7r/n — {m — l)Tr/n) -\- (m — l)Tx/n = tt — 2tx/ n — [m — l)Tx/n 

TT 

^ ALAjn+iA 2 ra+i = — (u—771—1). On the other hand due to the symmetry of the scaled polygon 
n 

(in red) ALAfn+iA 2 m+i = AA^+iMA' since AA^+iMN = AMAm+iA 2 m+i and ANMA' = 
AM Ajn+iL. AObMAm+i = 7r/2 {m — X)tx/ n as an exterior angle for l\MLA.m+\. 

Now, note that the scaled polygon towards the A^+i vertex (in red) reflected across the 
OfeL line segment will result in the scaled polygon towards the Am vertex (in blue); see also 
figure This is true because the ObL represents an axis of symmetry for the initial polygon 
and it contains point M which is a common point for both the scaled polygons towards the Am 
and Am+i vertices. Now we can deduce that these two scaled polygons will not intersect with 
each other if the vertices A' and A" stay together with Am+i on the same side with respect to 
the ObL axis. Then, AObMAm+i > AAm+iMA' and ALMAm+i > AAm+iMA" 
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Figure 4: Sketch of {n/m} star-polygon circumscribed in <Sb together with the scaled down 
polygon towards the Am+i vertex (in red) and its mirror image across the line (in blue) 
which is equivalent to the scaled down polygon towards the Am vertex. 


Am+l ^ AAm+lMA' 
■Kl2 - 1 - (m. — l)7r/n > tt — 2t: /n — {m — l)7r/n 
(4m — 4)/2n > (n — 4)/2n 
m > n/4 


ILMAm+i > lAm+iMA" 

■Kl2 — (m — l)7r/n > (m — l)7r/n 
1/2 > 2(m — l)/n 
n/4 -h 1 > m 


( 2 ) 


Thus, if the resulted fractal does not self-intersect, then n/4>m>n/4-|-l, which uniquely 
defines m except when n can be divided by 4 without residual. At the same time, if the strict 
inequalities of equations Q hold, this ensures that both of the vertices A' and A” stay on 
the right hand side of Oi,Lmie segment (see figure Q ^ their mirror images with respect to 
ObL will stay on the left hand side of O^L. Otherwise, if A' and A" were to cross the ObL 
segment becoming on its left hand side, the two polygons would intersect with each other and 
this intersection would be repeated everywhere since the resulted fractal is strictly self-similar. 
Finally, if A' or A" lies on OhL and the two scaled polygon have a common side MA' or MA”, 
then one of the equations (|^ must be with a sign for equality. 


P = 


Ml) 


2cos(^^^^)sin(^^f) 
P = Pi(n, m) for n = Av and m = v 

2cos(^^^^)sin(|^) 

2cos(^^^^)sin(f) 


Pi = 
Pi = 


Pi = 


Pi = 

Mi;) 


^/2cos((^) 


i2coAi)Mi) + V2Mi)Mi) 

Mi) 


Pi = 


Mi) 


■sin(^) 


P = P 2 (n, m) for n = Av and m = v -\-l 

Mi) 


P2 = 
P2 = 
P2 = 
P2 = 
P2 = 


2cos(ff)sin(h^) 

Mi) 

2 cos(f)sin(h^^) 


V^sin( 


(t; + l)7r 

Av 


Mi) 


V2sin(|f )cos(iL) + i2coAi) Mi) 
Mi) 


Mi) 


■sin(^) 


(3) 
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Therefore, n must be divided by 4 without residual because m G N. Thus, when n is divided 
by 4 without residual and the scaling ratio is P = P(n, m) from equations Q, then each of the 
scaled polygons has two common vertices with both of the adjacent scaled polygons. 

The case when n is divided by 4 without residual and m G [n/4, n/4+1] implies that equation 
0 produces two values for P. We will compute those values P(4n, v) and P(4n, n + 1), where 
n = Av for some n > 0,n G N; see equations Q. Equations Q clearly show that P(n, n/4) = 
P(n, n/4+1), when 4 divides n without residual, which ensures the unique definition of P(n, m) 
when the resulted fractals are non-self-intersecting. □ 

Corollary 1. The attractors of the IPS with n attracting points that are the vertices of {n,n/4} 
and {n,n/4+1} star-polygons and whose scaling ratios are P(n,n/4) and P(n,n/4+1), respec¬ 
tively, are identical and composed of star-polygons that have two vertices in common with each 
of the adjacent polygons. 

Remark 1. Note that theorem}^ is stated for the specific type of self-similar fractal sets ob¬ 
tained using IPS where the attracting points are the vertices of a {n/m} star-polygon and P 
is defined by n and m as given in equation 0- Thus, it does not exclude other kinds of 
non-self-intersecting star-polygonal fractal sets constructed in a different fashion. 


3.2 Fractal dimensions 


Theorem allow us to define P{n/m} as the IFS attractor produced from an initial equilateral 
{n/m}-polygon with n contracting maps that scale towards the n-vertices of the initial {n/m}- 
polygon with ratio P{n,m). We ensured that the obtained self-similar fractal set P{n/m} is 
non-self-intersecting, which allows us to compute its Hausdorff dimension dimHP{n/m} [2] by 
solving the following equation: 

N 



dimuFfnim} 


= 1 


(4) 


1=1 


where N indicates the amount of similarity maps Si (see Theorem and 0 < q < 1 are the 
scaling ratios for each similarity. 

Since P(2,1) = 1/2, dimHP{ 2 /i} = —ln{2)/ln{P{2, 1)) = 1 which means that the attractor 
of P{ 2 /i} is the initial line segment that connects both vertices or one can think about the 
Cantor set with scale ratio 1/2. 

For the Psgm in figure l^b) we have = 1 which lead to dimHP( 9 / 3 } = 

-H9)/ln(p/J)) « 1 6207686336597825. 

In the following figures ESIZI and 1^ are shown the attractors of P{n/m} where n = 
3,4,5,6,7,8,10,11,12,13,14,15,16,24 and m G [n/4,n/4 + 1]. The dimensions of the pre¬ 
sented fractals are computed in the examples that follows every figure. One can recognize well 
known fractals in the cases of n = 3,4, 5,6 but the other examples are not that famous due to 
the need of a special ratio in order to be constructed. 




Figure 5: IFS generated fractal sets made of points that lie in -F{ 3 /i}, -F'{ 4 / 2 }) -^{ 5 / 2 } and -^{ 6 / 2 } 
in subpanels (a), (b), (c) and (d), respectively. 


Example 1. IFe will compute the Hausdorff dimensions of the attraetors shown in figure\^- 


Figure 

Figure 

Figure 

Figure 


5(a) dimHF{ 3 /i} = -ln{3)/ln{P{3, 1)) 1.5849625007211563 

5(b) dimHF[ 4 / 2 } = —/n(4)//n(P(4, 2)) = 2 
5fc) dimHF{^/ 2 } = -/n(5)//n(P(5, 2)) 1.6722759381845547 

5(d) dimHF{ei 2 } = -/n(6)//n(P(6, 2)) 1.6309297535714573 
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Figure 6: IFS generated fractal sets made of points that lie in -F'{ 7 / 2 }) -^{ 8 / 3 }) -f{io/ 3 } and -Fpi/a} 
in subpanels (a), (b), (c) and (d), respectively. The IFS attractor of F^q/z} is shown in figure 

ib). 


Example 2. We will compute the Hausdorjf dimensions of the attractors shown in figure^ 


Figure 

Figure 

Figure 

Figure 


6 (a) dimHF{ 7 / 2 } = -ln{7)/ln{P{7, 2)) 1.6522616056918107 

6 (b) dimaF^s/s} = -ln{8)/ln{P{8,3)) ^ 1.6934291475411138 
6 (c) dimHF^w/ 3 } = -/n(10)//n(P(10, 3)) 1.5949906555938886 

6 (d) dimHF{ii/ 3 } = -/n(ll)//n(P(ll, 3)) 1.5911325154416658 
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Figure 7: IFS generated fractal sets made of points that lie in F{i 2 / 4 }, -^{ 13 / 4 }, -^{ 14 / 4 } and -^{ 15 / 4 }. 
in subpanels (a), (b), (c) and (d), respectively. 


Example 3. 

• Figure 

• Figure 

• Figure 

• Figure 


We will compute the Hausdorjf dimensions of the attractors shown in figure ^ 
'^a) dimHF{i 2 iA} = -/n(12)//n(P(12,4)) 1.598670034685813 

7(b) dimHF{i 3 / 4 y = -/n(13)//n(P(13,4)) 1.5653005271788485 

7(c) dimHF[u/A} = -/n(14)//n(P(14,4)) 1.5490615012592472 

7(d) dzmijP{i 5 / 4 } = -/n(15)//n(P(15,4)) 1.5430579163288531 
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Figure 8: IFS generated fractal sets made of points that lie in -F{i 6 / 5 } and -F{ 24 / 7 } in subpanels 
(a) and (b), respectively. 


Example 4. We will compute the Hausdorjf dimensions of the attractors shown in figure 


Figure 

Figure 


8 (a) dirriHFs^iQm = —/n(16)//n(P(16, 5)) 
lih) dimHF[ 24 / 7 } = -/n(24)//n(P(24, 7)) 


1.5434949184823248 

1.4772930562556852 


Also, the dimHF^n/m} for n G [17, 50] and m G [n/4, n/4 + 1] are as follows: 

{d7mijP{i7/5}., dzmj^F{i8/5},..., ~ (1.5238,1.5126,1.5071,1.5056,1.4924, 

1.4841,1.4794,1.4773,1.4677,1.4613,1.4573,1.4551,1.4478,1.4426,1.4391,1.437,1.4312,1.4269, 
1.4239,1.422,1.4172,1.4136,1.4109,1.4091,1.4051,1.402,1.3997,1.398,1.3946,1.3919,1.3898, 
1.3883,1.3853,1.3829} Finally, for n — le + 308, dimifP{ie+ 308 / 2 . 5 e+ 307 } ~ 1.001622. 


Theorem 3. As n goes to infinity, dimHF[n/m} approaches 1 


SlTiiTl /Ti] 

Proof. Let s = dimHFin/m\ then from lim P = lim ---r— -;- 

^ ^ n^oo n^oo 2cos{{m - 1)77/n)sin{m7r/ h) 


sin(w/n) 


= lim - 

n^oo 2 cOS (77/4) sm (77/4) 
ln(n) 

Thus, lim s = lim ——^ 

n^oo n^oo ln(? 7 / 77 ) 


= lim sin(K/n) = lim 77/n and nP^ = 1 we can deduce s. 


= 00 / 00 , hence lim s = lim 


d\n{n) 
dn 

n^oo n^oo ^ ln(n/7r) 

dn 


= 1 


□ 


As F{n/m} is inscribed in the same circle in which the initial {n/m}-polygon is inscribed, a 
corollary of Theorem is that as n —>■ 00 the F{n/m} is going to be arbitrary close to the circle 
in which the initial {^m}-polygon is inscribed. 


3.3 Exact drawing of the IFS iterations 

All the figures above were drawn by using a random walk generator that draws points which 
he in the IFS attractor EH- Another way of showing the attractor is by plotting large enough 
iteration (3th or 4th is usually enough) of the IFS where multiple scaled-down copies of the 
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Figure 9: IFS with initial {5/2} star-polygon where the second and the fourth iterations are 
shown in panels (a) and (b), respectively. 


initial polygon are imaged. In figure an example of this plotting approach is shown where 
in panel (b) the fourth iteration of the {5/2}-polygon looks like figure |^c) where the same 
attractor is produced by the random walk technique. In the next section we will use the latter 
technique more often for the sake of the clarity of the concepts presented. 


4 The centre of the circle as an additional point of at¬ 
traction 

If we add the centre of the circle as another attracting point that the random generator takes 
into account, then we can produce non-self-intersecting fractal sets that cover a great amount 
of the area that is bounded by the unit circle. This result is due to the fact that the centre point 
adds to the IFS attractor (the invariant set) one more scaled copy of the initial star-polygon, 
hence we need an additional contracting map which we will call Sc- Moreover, if the scaling 
factor of the central point is carefully computed, one can exploit a number of different features 
of the star-polygons. Here we will give a few introductory examples. 

Let us consider an initial polygon {3/1}, then P(3,1) = 1/2 and let us have a central 
map Sc with the same ratio 1/2 and rotation tt/S added to the set of maps {Si, S 2 , S 3 }. This 
will result in a triangular shape attractor with a Hausdorf dimension 7r/3] = 

—ln{4:)/ln{P{3, 1)) = /n(4)//n(2) = 2. Thus, the only difference from the attractor shown in 
figure l^b) will be the triangular shape. 

We will present the IFS of the {5/2} (see figure [lol Example with an attracting centre, 
where the similarity map corresponding to the centre point has the same scaling factor P(5, 2) as 
the maps that correspond to the vertices of the initial {5/2}-polygon. In figure [IT| another IFS 
is shown and its dimension is computed in Example]^ This fractal has centre-map that not only 
scales, but also rotates the initial polygon at angle 7r/5 while keeping the ratio P(5, 2). Let us 
also see the lES of the {6/2} with an attracting centre, where the similarity map corresponding 
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to the centre point has a scaling factor P(6, 2) and does not imply rotation, presented in figure 
T^and Example!^ 

From the dimensions computed in section |3.2| we can deduce that as n grows the scaling 
ratio P(n,m), where m G [n/4, n/4+ 1] monotonically decreases and the resulting fractals 
shrink in dimension. Thus, if we want to increase n, but keep the attractors with a reasonably 
high dimension, we can no longer use the same ratio for the centre scaling map as in the cases 
for n = 5 and n = 6. Therefore, we will define different rules for the scaling ratio of the centre 
map Sc for any n, depending if it is odd or even and if Sc includes any rotation such as yr/n or 
it does not. 



Figure 10: IFS with initial {5/2} star-polygon where the central map Sc has no rotation and 
uses the ratio P(5, 2). The first, the second and the fourth iterations are shown in panels (a), 
(b) and (c) respectively. 100000 points that he on the attractor of the IFS are shown in panel 
(d). 
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Example 5. We will compute the Hausdorjf dimensions of the attractor shown in figure\ld(d): 
dimHF{^l 2 }[L^A = -/n(6)//n(P(5,2)) 1.8617 



Figure 11: IFS with initial {5/2} star-polygon where the central map Sc has rotation 7r/5 and 
uses the ratio P(5, 2). The first, the second and the fourth iterations are shown in panels (a), 
(b) and (c) respectively. 100000 points that he on the attractor of the IFS are shown in panel 
(d). 


Example 6. We will compute the Hausdorff dimensions of the attractor shown in figure\Tl^d): 
dimHF^5/2}[L^,Tv/5] = 0] = —/n(6)//n(P(5, 2)) ~ 1.8617 
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Figure 12: IFS with initial {6/2} star-polygon where the central map Sc has no rotation and 
uses the ratio P(6, 2). The first, the second and the fourth iterations are shown in panels (a), 
(b) and (c) respectively. 100000 points that lie on the attractor of the IFS are shown in panel 
(d). 


Example 7. We will compute the Hausdorff dimensions of the attractor shown in figure\T§(d): 
dimHF{e/2}[L^,0] = -ln{7)/ln{P{6,2)) ^ 1.7712 
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Figure 13: Sketch of a {7/2} polygon IFS with a centre map after the first iteration. Shown is 
the way of computing of the ratio for the centre map. 


In order to clearly show how the ratios of Sc are deduced, in figure [T^ is sketched part of 
the IFS of {7/2}-polygon when it is iterated only once. There are two centre polygons, one 
is rotated at angle n/7 (the one with a vertex at M^) and another rotated at an angle to be 
computed later (with vertex L^). Thus, the contraction map that scales the original {7/2}- 
polygon towards the point O and images it in one of these two polygons includes rotation as 
well. Using the sketch in figure [I^ we will show that some of the points M* and L* will be used 
as vertices of a central polygon that corresponds to a central map Sc, which can be used for 
non-self-intersecting polygonal IFS. 

Let us construct the points M^, M^, M^,..., lying on the line that crosses the line- 

segment 0^0 at angle T^jl. This line is intersected by the line segments that start from 
the vertices L},L^ and and are perpendicular to the line-segment 0^0 resulting in the 
points M^,M^ and M^. Now we will define six different ratios Pc for the centre map Sc'. 
OM^/OA, OM^/OA, OM^jOA and OL^/OA, OV'jOA, OL^/OA. All of them are de¬ 
fined by the angles — It^/u where in general I — 2z -|- 1 for i > 0,z G Z when 

the polygons are odd-sided and I = 2i when the polygons are even-sided. For I = 1,3,5, 
ON = 0 A- O^A - cos(/7r/7) and OM^ = ON/ cos(7r/7), thus OM^/OA = 1/ cos(7r/7) - 

{O^A/OA){ 1 /cos{ 7 t/ 7 ) — cos(/7r/7)/cos(7r/7)). This equation also holds for I = 7, where 
A = LS and if apply equation Q, so that O^A/OA = P{7, 2), then OM^/OA = 1/ cos(7r/7) — 
P(7,2)/cos(7r/7) — P(7, 2) cos(/7r/7)/cos(7r/7). Also, if generalised for an arbitrary initial 
{n/m}-polygon, it leads to the following equation: 


OM^ 1 — P{n,m) — P{n,m) cos{l 7 r/n) 

OA cos(7r/n) 

AO^OM^ = Ti/n 

0 < I < n,l = 2 i if n-even, / = 2z -|- 1 if n-odd, i > 0 ,i E Z 
Another ratio that may be used for the map Sc is OL//OA. Here sin(7r/7). 
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hence, tan(ZO^OL^) 


- 


sin(7r/7) 


OA-O^A-O^L^ cos(7r/7) 
sin(/7r/7) 


which for an arbitrary I will become 


tan(ZO^OL^) = — -——- ^, --—7—. Now let us take into account that = O^A 

^ ' OA-O^A-O^Ucos{l7r/7) 


and O^A/OA = P(7,2), hence tan{AO^OL^) = ^ 2)’-^P(7 


o^ n ~FA- Finally, for 
2 ) cos(i7r/7) 

an arbitrary initial {n/m}-polygon we can deduce the angle AO^OL} and from {OVY — 
{NVy + {ONY we can also deduce the ratio OL)/OA as follows: 


OL' 


, m){P{n, m) — 1)(1 + cos(/7r/n)) + 1 


( 6 ) 


/ 7 X ( P(n,m) sin(/7r/n) 

7 (n, m, 1 ) = AO OL^ = arctan --—-r-—-^- 

\ 1 — -T (n, m) — P[n, m) cos[lw/n) 

0 < I < n,l = 2i if n-even, / = 2i + 1 if n-odd, z > 0, z G Z 


Now we can look back at figures nnni and [ 1 ^ and understand how they are constructed. 
In example]^ figure 10, the contraction ratio of Sc is OL^/OA, and the attractor of the IFS 
is denoted as F{5/2}[7/\ 0], where 0 indicates the angle of rotation that Sc has. In this case 
of {5/2}-polygon OM^/OA = OL^/OA, so it does not matter if or is used for the 
notation. In the other examples where OV'/OA = OM^ jOA^ again iJ will be used as notation. 
In example]^ figure [T^ the contraction ratio of Sc is again OL}/OA, but here we have rotation 
at angle 7r/5, thus the attractor of the IFS is denoted as F{ 5 / 2 }[F^, 7r/5]. And finally in example 
figure \A2\ the contraction ratio of Sc is OlAjOA and the attractor of the IFS is denoted 
F^{6/2}[A°,0]. 


4.1 Even n 


In this subsection we will take a close look at the IFS that originates from even sided star- 
polygons. Firstly, we should note that the same way as F{ 6 / 2 }[A°, 0], for any even n, F{ 2 i/m}[A°, 0] 
will always be a non-self-intersecting attractor if m G [n/4, n/4 -|- 1] and z G N. Now, the first 
example has {6/2} as initial polygon, and Sc has scaling ratio OLA/OA and rotation tt/G. 
The resulting fractal can be seen in figure where from the random generated attractor, see 
panel (d), we can expect the exact dimension of 2. Indeed, this is analytically proven in the 
computations of example Similarly, the constructions and the attractors of F'{ 8 / 2 }[A^ |] and 


F{ 8 / 3 }[T^, |] are shown in figures 15 and 16[ They have equal dimension computed in example 
1^ Another pair of attractors that have central map and originate from {8/2}-star polygon 


are the F{8/2}[T°,0] and F{8/3}[T°,0] shown in figures 17 and 18 They have equal dimension 


computed in example The last four examples of attractors clearly show that the scaling 
ratio and the number of the vertices are the parameters that define the attractor of the IFS. 
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Figure 14: IFS with initial {6/2} star-polygon where Sc has scaling ratio OL^/OA and rotation 
tt/6 . The hrst, the second and the fourth iterations are shown in panels (a), (b) and (c) 
respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 


Example 8. We will compute the Hausdorff dimensions of the attraetor shown in figure\If^d): 

+ y/ 2 P( 6 , 2 )(P( 6 , 2) - 1)(1 + cos(2ir/6)) + = i 

Hence, QiS + y — 1 = Q and Vi 2 = 1/3; —1/2, therefore as y > 0 

^ 1/3 ^ |] = 2 . 
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Figure 15: IFS with initial {8/2} star-polygon where Sc has scaling ratio OL^/OA and rotation 
tt/8 . The hrst, the second and the fourth iterations are shown in panels (a), (b) and (c) 
respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 


Example 9. The Hausdorjf dimensions of the attractor shown in figure\r^d) is 
dimHF{s/ 2 }[L‘^, |] ~ 1.9799 
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Figure 16: IFS with initial {8/3} star-polygon where Sc has scaling ratio OL^/OA and rotation 
tt/8 . The hrst, the second and the fourth iterations are shown in panels (a), (b) and (c) 
respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 
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Figure 17: IFS with initial {8/2} star-polygon where Sc has scaling ratio OL^/OA and no 
rotation. The first, the second and the fourth iterations are shown in panels (a), (b) and (c) 
respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 


Example 10. The Hausdorff dimensions of the attractor shown in figure\P^d) is 
dimHF{s/2}[L^,0] ~ 1.8678 
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Figure 18: IFS with initial {8/3} star-polygon where Sc has scaling ratio OL^/OA and no 
rotation. The first, the second and the fourth iterations are shown in panels (a), (b) and (c) 
respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 


4.2 Odd n 


In this subsection we will take a close look at the IFS that originates from odd-sided star- 
polygons. The first example has {7/2} as an initial polygon and Sc has a scaling ratio OM^/OA 
and rotation tt/T; see Eqs. Q. The resulting fractal can be seen in figure 19, where the 
first, the second and the fourth iterations are in panels (a), (b) and (c), while the randomly 
generated attractor is in panel (d). The Hausdorf dimension of the attractor ,tt/7] 

is computed in example [TT| to be 1.8773. Another fractal that originates from a {7/2}-polygon 
is shown in figur e Here the scaling ratio is OL^/OA and the angle of rotation is computed 
using equations (1^, which leads to polygons that meet at their vertices. The dimension of 
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7 ( 7 , 2,1)] is computed in example 12 to be 1.8564. 


The technique that uses the ratio OL^/OA and the angle from equations ^ can also be 
used for producing non-intersecting self-similar fractals for any odd n. Therefore, for any n we 
can generate a n-flake which will be either F{n/m}[L^, 0 ] if n is even or 7 ( 72 , m, 1 )] if 

n is odd. 


4.3 The dimension of 00 -flake 

For any n, m G [n/4, n/4-|-1] and z G N, by construction ii n = 2i then F{n/m}[L^, 0] is non-self- 
intersecting, and by construction if n = 2i-|-1 then F^n/m} [L^ , li'n, m, 1)] is non-self-intersecting. 
Therefore, equations Q, and Q can be used for the corresponding dimHF^n/m}[L ^and 
dimHF^n/m}[L^ l)]to be obtained. 

Theorem 4. As n goes to infinity, dimHF{n/m}[LO,o] o,nd dimHF^n/m}[L^approach 2. 
Proof. Both dimensions can be deduced from the equation 


nP{n, my F \/2P{n, m){P{n,m) — 1){1 F cos{h^/n)) + 1 = 1, 


where s denotes dimHF^n/m}[L°fl] or dimHFi^n/m}[L^,Fn,m,i)\- The latter equation can be modified 
to 

p. = i _ 7^ , m){P{n^ m) — 1)(1 + cos(/7r/n)) + 
n n 


from where: 


1 y/^ , m){P{n, m) — 1)(1 -|- cos(/7r/n)) -|- 1 


s = lim 

n—^-oo 


n 


n 


ln(P) 

1 y/jPG , m){P{n, m) — 1) F 1 


= lim 

n—^-oo 


n 


n 


ln(P) 


In (1 - (1^ 


= lim 

n—^-oo 


n 


n 


1 (l-2n/ny 


HP) 

Let us substitute v = t; jn, hence 


= lim 

n^oG 


n 


n 


ln(7r/n) 


lnf-(l-(l-20®^ 

s = lim- - -——-^ = lim , ..... 

i^->o ln(i/) ln(z/) I'-i-o ln(z/) 


In ( - 

TT 


-|- lim 


ln(l - (1 - 20") 


-oo 2si/(l-20"“^ ,. 2si/(l-(s-l)2i/ + 0(2)) 

“ ^ ^ ^ ^ ^ 1 - (1 - 2uY “ ^ ^ 1 - {1 - s2u F 0{2)) 


, ,. 2 siyFO(2) ^ 

= 1 F hm- , = 2 

2 siy -j- 0(2) 


□ 
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Figure 19: IFS with initial {7/2} star-polygon where Sc has scaling ratio OM^/OA and rotation 
tt/7. The first, the second and the fourth iterations are shown in panels (a), (b) and (c) 
respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 


Example 11. The Hausdorff dimensions of the attractor shown in fiqureUWd) is 
dimHF{y2 }[M\^] ^ 1.8773 
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Figure 20: IFS with initial {7/2} star-polygon where Sc has scaling ratio OL^/OA and rotation 
7(7, 2,1) radians. The first, the second and the fourth iterations are shown in panels (a), (b) 
and (c) respectively. 100000 points that lie on the attractor of the IFS are shown in panel (d). 


Example 12. The Hausdorff dimensions of the attractor shown in figure^dfd) is 
dimfi-F{7/2}[Li,7(7,2,i)] ~ 1.8564 

5 Special cases and equivalent IFS attractors 

As in the previous sections, here we assume n > 2, n e Z, m e Z, 0 < m < n, and z G Z, z > 0. 
Let us now review the developed notation and how many of the well known fractals can be 
associated with it. Firstly, we can say that the Cantor set results as an IFS attractor if rz = 2 
and P = 1/3. As we saw in figure [^a), the Sierpinski Triangle comes out when n = 3 and 
P = P(3,1) = 1/2 or P{3/i} and the Sierpinski Hexagon when n = 6 and P = P(6, 2) = 1/3 or 
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Figure 21: In subpanels (a) and (b) we can see the third iteration of the IFS that originate 
from {6/3} and {8/4}, respectively. Both Sc have scaling ratio OL^/OA, where in panel (a) 
the rotation is tt/G and in (b) it is tt/S. In subpanels (c) and (d) we can see the third iteration 
of the IFS that originate from {7/3.5} and {9/4.5}, respectively. Both Sc have scaling ratio 
OM^/OA, where in panel (c) the rotation is tt/T and in (d) it is tt/Q. 
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F{ 6 / 2 }- The Greek Cross fractal appears when n = 4 and P = P(4, 2) = 1/2 or F{ 4 / 2 }, while for 
n = 4 and P < lf2 the invariant set for a Horseshoe map is produced. The Sierpinski Pentagon 
appears for n = 5 and P = P(5, 2) = 1/(1 + golden ratio) or P{ 5 / 2 }- When the centre map Sc is 
taken into account, the Vicsec fractal can be produced when n = 4 and P = Pc = 1/3. Also, the 
Pentaflake and the Hexaflake are shown in figures 11 and 12 as F{5/2}[T^,0] and P{ 6 / 2 }[T°, 0], 
respectively. 

However, the attractors of the special cases mentioned above do not originate from an unique 
star-polygons because if we generate random points on the attractor of the IPS, the image is 
defined by the maps {A*}; see theorem Therefore, if we do not alter P and n, the IFS- 
attractors with any initial {n, m}-polygon will be the same. Thus, we can assume that every 
attractor that originates from a {2i,m}-polygon is equivalent to the attractor of the IPS that 
originate from the {2i,i}-polygon, where the ratios and the rotations of the maps {Si, ...,S 2 i, Ac} 
are kept the same as the ones used for the IPS of the {2i,m}-polygon. However, when the exact 
polygons are plotted, as in section]^ due to the impossibility for infinite iterations to be realised, 
the integer m also plays its role. Two such examples are shown in figure where in panels (a) 
and (b) the third iterations of the IPS of F| 6 / 3 >[T^, |] and F{ 8 / 4 }[L^, |] are realised respectively. 


One can see the difference with figures and 

F{s/3}[L‘^, |] are shown. 
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where F{ 6 / 2 }[T^f], F{ 8 / 2 }[T^|] and 


In the case of odd n = 2f + 1 we do not have a star-polygon composed of n line-segments 
that cross each other at the centre and are inscribed in all the {2i-|-l,m}-polygons in the same 
way as the {2i,i}-polygon is inscribed in any {2i,m}-polygon. Therefore, we will construct 
such polygons as the star of 2z -|- 1 line segments that start from the centre of a 2i -|- 1 regular 
polygon and end up at its vertices. Let us denote this figure as {2i -|- 1, f -|- l/2}-star polygon. 
As the {2i -|- l,f -|- l/2}-polygon is inscribed in all the {2i -\-1 , m}-polygons we can generalise 
that any attractor that originates from a {n,m}-polygon is equivalent to the attractor of the 
IPS that originates from the {n,n/2}-polygon where the ratios and the rotations of the maps 
{Ai,..., Sn, Ac} are kept the same as the ones used for the IPS of the {n,m}-polygon. Two such 
examples are shown in figure where in panels (c) and (d) the third iterations of the IPS of 
-A{ 7 / 3 . 5 }[M^, |] and F{ 9 / 4 . 5 }[M^, |] are realised. One can see the difference with figures [l^ where 
T{ 7 / 3 }.[M\ y] was shown. 


Example 13. If the fractal figure \2l\/d) is infinitely iterated the attraetor will have Hausdorff 
dimension: dimHF{ 9 / 4 . 5 }[M^, |] ~ 1.8879 


Unlike the cases of n = 3,5,7 and 9, where the is a non-self intersecting 

fractal with the copies of the initial {n/|}-polygon osculating with each other, for n = 2f -|- 1 
when n > 11, the scaled copies stop osculating and with the increase of the iterations they 

do not fill up the space in the most effective way. The same effect appears when we take 

for n — 2i when n > 10. Thus, the problem of finding scaling ratio for the Ac for 

every n > 10 where the rotation of Ac is equal to 7r/n stays as an open problem. This is an 

important question due to the fact that the fractals that originate from an n-gon with rotation 
of Ac equal to ir/n could have dimensions very close or equal to 2; see examples and 


Conclusion 

The present paper develops a universal technique that allows any star-polygon to be used for the 
construction of non-self intersecting fractal (Sierpinski n-gon or n-fiake) by using IPS through 
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random walk or through an exact scaling. Along the proposed scaling ratios, the Matlab code 
for IFS random walk fractal generation is provided so that someone interested in studying the 
geometry of this class of fractals could use it. Important dimensions are computed, namely 
the dimension of the Sierpinki oo-gon is proved to be 1, the dimension of the oo-flake to be 2, 
and the dimension of F{6, 2}[L^, 7r/6] to be 2 as well. It is also shown, that by using random 
walk IFS generator, identical attractors may result from different initial star-polygons. The 
proposed study can be extended if rotations are applied not only to the Sc map, but to the Si 
maps as well. However, this is still an ongoing research in development. 

The techniques for construction and the provided ratios needed for the dimensions of the 
presented class of fractals is important not only for mathematicians, but also for engineers and 
other scientists who may be interested in fractal-shaped devices or who study the fractal shapes 
of nature. With the advanced precision of the fabrication technology polyflakes may become 
an important design for devices such as antennas or chemical mixers for fuel cells, batteries, 
etc. Also, fractal shapes are applicable in any kind of wave absorbers where the wave could 
be a sound wave, electromagnetic signal, light, caused by turbulent flow, etc. In other words, 
fractal designs are going to be part of the future physical devices at all scales and, hence, the 
research focused on fractal-shaped figures is important for many innovation processes. 
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Appendix (Matlab random generator) 


% name: Vassil Tzanov 


% function ’frac’t^ 


hat takes the matrix C, 


%C 


A1 bl 
A2 62 


\An bn J 

% Ai are matrices 2x2 ,bi are vectors 2x1, which are the i-th linear function in the ”IFS” 
% k defines the amount of points that we want to be plotted 
% on the attractor defined by C; the function ’frac’ can plot fractals 
% with central polygon rotated at angle rot; the central map must be defined by 
% the last two rows of (7; if there is no central map, rot has to be defined as 0 
function Fractal = frac(C,k,rot) 


% verification of C 
D=size(C); 
n=D(l)/2; 

if D(2)~=3 I floor(n)~=n 
Fractal = ’Bad input size’; 
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else 


% computation of the determinants of the matrices A 

for i=l:n 

dets(i)=abs(det(C((2*i-l);(2*i),l;2))); 
end 

dets=dets’; 

% if some determinant is zero we have to add a small value, because 
% we do not like this function to be executed with zero probability 
dets=max(dets, max(dets)/(25*n)); 

% the determinant are divided on their sum so we derive a probability vector 
dets = dets/sum(dets); 

% the vector prob is dehned as 

% prob = [0, dets(l), dets(l)+dets(2), . . dets(l)+...+dets(n-l)] 

% then ” sum (prob < rand)” gives a random integer between 1 and n 
% that we use in order to randomly choose one of the n-th functions 
de=dets(l); 
prob(l)=0; 
for i=2:n 
prob(i)=de; 
de=de+dets(i); 
end 

% computations of the points’ coordinates and hlling the matrix Fractal 
x=[0;0] ; 

Fractal=zeros(2,k+20); 
for j=l:(k+20) 

int=sum(prob<rand); 

% one of the functions is randomly chosen 
i=int; 

% matrix computation of t * x + (1 — t) * y, where x is 
% the current point and y is the point towards which x is attracted 
% if z corresponds to the last row, 

% a rotation at rot radians around [0,0] must be executed 

if i==n 

x=[cos(rot),-sin(rot);sin(rot),cos(rot)]*x; 

x=C((2*i-l): (2*i)4:2)*x+(diag([l,l] ,0)- 

C((2*i-1):(2*i),l:2))*C((2*i-l):(2*i),3); 

Fractal(:,j)=x; 
else 

x=C((2*i-l):(2*i),l:2)*x+(diag([l,l], 0 )- 

C((2*i-1):(2*i),l:2))*C((2*i-l):(2*i),3); 

Fractal(:,j)=x; 
end 
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end 


% the points of the matrix Fractal are plotted 

Fractal=Fractal(:,21:(k+20)); 

plot(Fractal(1,;).Fractal(2,:),’k. \ ’MarkerSize’,1) 
axis Coequal’) 
end 
end 
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